o 

oo 



> 

0^ 



O 



B 

> 

X 



Relaxation of Excited States in Nonlinear 
Schrodinger Equations 

Tai-Peng Tsai* Horng-Tzer Yau^ 
Courant Institute, New York University 



O : August 27, 2001 

O 
(N 



Abstract 



We consider a nonlinear Schrodinger equation in M with a bounded local 
potential. The linear Hamiltonian is assumed to have two bound states with the 
eigenvalues satisfying some resonance condition. Suppose that the initial data 
I is small and is near some nonlinear excited state. We give a sufficient condition 

■ on the initial data so that the solution to the nonlinear Schrodinger equation 

approaches to certain nonlinear ground state as the time tends to infinity. 



i!B ■ 1 Introduction 



Consider the nonlinear Schrodinger equation 

tdt^ = (-A + V)^ + X\^\^ij, ^{t = 0)=^o (1.1) 

where is a smooth localized potential, A is an order 1 parameter and ip = il){t^x) : 
M X M'^ — > C is a wave function. The goal of this paper is to study the asymptotic 
dynamics of the solution for initial data ip^ near some nonlinear excited state. 
Recall that for any solution il){t) G H^{M?) the L^-norm and the Hamiltonian 

=Jlm' + lv\^\' + \xm'dx , (1.2) 

are constants for all t. The global well-posedness for small solutions in if^(]R^) can 
be proved using these conserved quantities and a continuity argument. 
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We assume that the hnear Hamiltonian Hq := —A + V has two simple eigenvalues 
Co < ei < with normalized eigen-functions (po, 0i. We further assume that 



eo<2ei. (1.3) 



The nonlinear bound states to the Schrodinger equation (|LT|) are solutions to the 
equation 

{-A + V)Q + X\Q\^Q = EQ . (1.4) 

They are critical points to the Hamiltonian defined in ( |1.2|) subject to the con- 
straint that the L^-norm of is fixed. For any bound state Q = Qe, 'ipit) = Qe~^^^ 
is a solution to the nonlinear Schrodinger equation. 

We may obtain two families of such bound states by standard bifurcation theory, 
corresponding to the two eigenvalues of the linear Hamiltonian. For any E sufficiently 
close to Co so that E — cq and A have the same sign, there is a unique positive solution 
Q = Qe ( |1.4| ) which decays exponentially as x — > cxo. See Lemma pri| . We call this 
family the nonlinear ground states and we refer to it as {Qe}e- Similarly, there is a 
nonlinear excited state family {Qi^e}e- We will abbreviate them as Q and Qi. From 
Lemma 01 we also have \\Qe\\ ^ \E — CqI^^"^ and ~ l-^' ~ ei| 
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It is well-known that the family of nonlinear ground states is stable in the sense 
that if 



inf ||^(t)-Q^e 



ie\ 



Il2 



is small for t = 0, it remains so for all t, see H]. Let ||-||r2 denote a local norm 

loc 

( a precise choice will be made later on). One expects that this difference actually 
approaches zero in local norm, i.e.. 

Mm miU{t) -Qec'^L, =0. (1.5) 

t^OO &,E " "-'^loc 

If —A + V has only one bound state, it is proved in [|l^ ^ that the evolution will 
eventually settle down to some ground state Qeo^ with Eoo close to E. Suppose now 
that —A + V has multiple bound states, say, two bound states: a ground state (po 



with eigenvalue Cq and an excited state 0i with eigenvalue ei. It is proved in |12] that 
the evolution with initial data tpo near some Qe will eventually settle down to some 
ground state Qeoc with Eoo close to E. See also [|I| for the one dimensional case and 



TT| for nonlinear Klain-Gorden equations. 
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Denote by Lf. the weighted spaces (r may be positive or negative) 

LliW") = {0eL2(R3) : (x)"0gL2(R3)} . (1.6) 

The space for initial data in is 

Y = H\R^) n L^^ (M^) , ro > 3 . (1.7) 

We shall use Lf^^ to denote L'^ro- The parameter ro > 3 is fixed and we can choose, 
say, ro = 4 for the rest of this paper. We now state the assumptions in [|T^ on the 
potential V. 

Assumption AO: —A + V acting on L^(R^) has 2 simple eigenvalues eo < ei < 0, 
with normalized eigenvectors (po and (pi. 

Assumption Al: Resonance condition. Let eoi = ei — eo be the spectral gap of the 
ground state. We assume that 2eoi > |eo| so that 2eoi is in the continuum spectrum 



of Hi. Let 



70 := lim Im Lo^?, tTT \ : Pc^'Vo^? 

0-^0+ \ Hq + eo — 2ei — cri 



Since the expression is quadratic, we have 70 > 0. We assume, for some sq < Cn^ 
small enough, 

inf hm Im f 0o0?, tt^ : Pc^'Vo^?) > ^70 > . (1.9) 

We shall use Oi to replace ai and the limit limo-^o+ later on. 

Assumption A2: For \Q\ sufficiently small, the bottom of the continuous spectrum 
to —A + V + \Qe, 0, is not a generalized eigenvalue, i.e., not a resonance. Also, we 
assume that V satisfies the assumption in [|1^ so that the W'^'^ estimates k < 2 for 



the wave operator Wh = lim^^oo e'^^e***-^"'"'^-' hold for k < 2, i.e., there is a small 
a > such that, 

|V°\/(x)| < for|a|<2. 

Also, the functions (x ■ V)^V, for k = 0, 1,2,3, are —A bounded with a — A-bound 
< 1: 

||(x- V)'^y0||2 <ao||-A0||2 + C||0||2, (To < 1, A; = 0,1, 2, 3 . 

Assumption A2 contains some standard conditions to assure that most tools in 
linear Schrodinger operators apply. These conditions are certainly not optimal. The 
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main assumption in A0-A2 is the condition 2eoi > |eo| in assumption Al. The rest of 
assumption Al are just generic assumptions. This condition states that the excited 
state energy is closer to the continuum spectrum than to the ground state energy. It 
guarantees that twice the excited state energy of Hi (which one obtains from taking 
the square of the excited state component) becomes a resonance in the continuum 
spectrum (of Hi). This resonance produces the main relaxation mechanism. If this 
condition fails, the resonance occurs in higher order terms and a proof of relaxation 
will be much more complicated. Also, the rate of decay will be different. 

The main result in |jl2| concerning the relaxation of the ground states can be 



summarized in the following theorem. 

Theorem A Suppose that suitable assumptions on V hold. Then there are small uni- 
versal constants eo,nQ > such that, if the initial data ipo satisfies W^po ~ ""-e^^^^oHy ^ 
eln^ for some n < uq and some Qq G M, then there exists an Eoo and a function G(t) 
such that WQeocWy — n = 0{eln), &{t) = —E^t + O(logt) and 

\\m-QE^e^''^'^\^. <c(l+t)-^/^ (1.10) 

loc 

This theorem settles the question of asymptotic profile near ground states. Sup- 
pose that the initial data ipQ is now near some nonlinear excited state. From the 
physical ground, we expect that ipt will eventually decay to some ground state un- 
less the initial data ipQ is exactly a nonlinear excited state. We call this the "strong 
relaxation property". For comparison, we define a weaker property, the "generic re- 
laxation property", as follows. Denote the space of initial data by X. Let Xi (Xq 
resp.) be the subspace of initial data such that the asymptotic profiles are given by 
some nonlinear excited (ground resp.) states. We shall say that the dynamics satisfy 
the generic relaxation property if Xi has "measure zero". This concept depends on 
a notion of measure which should be specified in each context. 

With this definition, the strong relaxation property means that Xi is exactly the 
space of excited states. In particular, Xi is finite dimensional. We first note that the 
strong relaxation property is false. For any nonlinear excited state Qi, define Xi^q^ 



to be the set of initial data converging to Qi asymptotically. It is proved in ||T3| that 
for any given nonlinear excited state Qi, Xi^q^ contains a finite co- dimensional set. 
Thus our goal is to establish some weaker statement such as the generic relaxation 
property. This is the first step toward a classification of asymptotic dynamics of the 
nonlinear Schrdinger equation. 
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In this paper, we shall prove that for any excited state, there is a small neighbor- 
hood N so that 

lATnXil < cU4^\N'nXo\ 

This estimates states that the ratio between Xi and X around excited states are 
bounded by the mass of the initial wave function HV^oP- (Since we have not given a 
measure on the space of initial data, this statement is not well-defined and a precise 
statement will be given later on.) 

In order to state the main result, we first decompose the wave function using the 
eigenspaces of the Hamiltonian Hq as 

V' = x(^o + y</'i + i, e=Pc''''V'- (1-11) 

For initial data near excited states, this decomposition contains an error of order 
and it is difficult to read from (|1.11| ) whether the wave function is exactly an excited 
state. Thus we shall use the decomposition 

^ = x0o + gi(2/)+e- (1-12) 

where 

y = y_, x = x-(</>o,Qi(2/)), i = i- V^y) . (1.13) 

Here we have used the convention that 

Qi{y) ■■= Qiim)e^ , m=\yl me'® = y 

We shall prove that for with sufficiently small Y norm ( [1.7D , such a decomposition 
exists and is unique in section 2. Thus we assume that ipo = Xo0o + Qiivo) + is 
sufficiently small in Y. Let and be the same small constants given in Theorem 
A. By choosing a smaller no, we may assume uq < Eq/A. We assume the initial data 
satisfies that 

llV'olly = n, <n<no, \yo\ > \n, 

(1.14) 

\xo\ >2ne |xo| > e^^n'^ ||^o||y , 

where £2 > is a small universal constant to be fixed later in the proof. These 
conditions can be interpreted as follows: The excited state component, yo, should 
account for at least half the mass of the initial data. Under this condition, if the 
ground state component, Xq is not too small compared with the continuum component 
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^, then the dynamics relaxes to some ground state. The condition |xo| > 2ne~^/^ is 
a very mild assumption to make sure that Xq is not incredibly small. The following 
constant will be used in many places in this paper. Define 

e := min (log(2n/xo))-'/'} . (1.15) 

Since n <nQ< e^/A and |a;o| > 2ne^^^"', we have n < e^. 

Theorem 1.1 Suppose the assumptions on V given above hold. Let ip{t,x) be a 
solution of (pTTP with the initial data ipo satisfying (|1.14|) . Let 



n, = {\xo\' + llvol'Y^' - n . (1.16) 

Then, there exists an and a function Q{t) such that HQsoolly ~ ""-i = 0{en), 
0(t) = -E^t + O(logt) and 



C,{1 + t)~V2 < ll^(^) _ g^^e^®^*^ IL. <C,{1 + t)-"^ . (1.17) 

loc 

for some constants Ci and C2. 

It is instructive to compare our result with the linear stability analysis of 
0] and 0. In our setup the main result in states that the linearized 
operator around a nonlinear excited state is structural stable if cq < 2 ei and unstable 
if Co > 2ei. Hence the excited states considered in this article is unstable and 
is expected to decay under generic perturbations. The instability of the excited 
state contains in Theorem 1.1 is thus consistent with the linear analysis. Notice 
that Theorem 1.1 tracks the dynamics for all time including time regime when the 
dynamics are far away from the excited states. Furthermore, for all initial data 
considered in Theorem 1.1, the relaxation rate to the asymptotic ground state is exactly 
of order a rate very different from the standard linear Schrodinger equations. 

In view of the linear analysis, the existence |TB[ of (nonlinear) stable directions for 



excited states is a more surprising result. For the linear stable case, i.e., cq > 2 ei, the 
only rigorous result is the existence [|l^ of (nonlinear) stable directions in this case. 



Although the linear analysis states that all directions are linearly stable, on physics 
ground we still expect excited states remains generically unstable. 

We now explain the main idea of the proof for Theorem 1.1. The relaxation 
mechanism can be divided into three time regimes: 
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1. The initial layer: The component of the wave function in the continuum 
spectrum direction gradually disperses away; the components in the bound states 
directions do not change much. 

2. The transition regime: Transition from the excited state to the ground 
state takes place in this interval. The component along the ground state grows in 
this regime; that along the excited state is slightly more complicated. We can further 
divide this time regime into two intervals. In part (i), the component along the excited 
state does not change much. In part (ii), it decreases steadily and eventually becomes 
smaller than the component along the ground state. 

3. Stabilization: The ground state dominates and is stable. Both the excited 
states and dispersive part gradually decay. 

In different time regimes, the dominant terms are different and we have to linearize 
the dynamics according to the dominant terms. In the first time region, ^/'(t) is near 
an excited state, and it is best to use operator linearized around the excited state. 
In the third time regimes, iplt) is near a ground states, and it is best to use operator 
linearized around a ground state. For the transition regimes, the dynamics are far 
away from both excited and ground states and we can only use the linear Hamiltonian 
Ho. 

Besides technical problems associated with changing coordinate systems in dif- 
ferent time intervals, there is an intrinsic difficulty related to the time reversibility 
of the Schrodinger equation. Imagine that we are now ready to show that our dy- 
namics is in the third time regime and will stabilize around some nonlinear ground 
state. If we take the wave function ipt at this time and time reverse the dynamics, 
then the dynamics will drive this wave function back to the initial state near some 
excited state. The time reversed state ipt and the wave function ipt itself will satisfy 
the same estimates in the usual Sobolev or Lp senses. However, their dynamics are 
completely different: one stabilizes to a ground state; the other back to near an ex- 
cited state. This suggests that ipt carries information concerning the time direction 
and this information will not show if we measure it by the usual estimates. 

This time reversal difficulty manifest itself in the technical proofs as follows. We 
shall see that, when the third time regime begin, the dispersive part is not well- 
localized and its L^-norm can be larger than that of the bound states — both violate 



conditions for approaching to ground states in |T^. To resolve this issue, we need 
to extract information which are time-direction sensitive so that even though the 
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disperive part may be large, it is irrelevant since it is "out-going". Though the 
concept of "out-going" wave is known for linear Schrodinger equations, it is difficult 
to implement it for nonlinear Schrodinger equations. Our strategy is to identify the 
main terms of the dispersive part and calculate them explicitly. These terms carry 
sufficient information concerning the time direction. The rest are error terms and we 
can use various Sobolev or Lp estimates. 

Resonance induced decay and growth 

To illustrate the mechanism of resonance induced decay and growth, we consider 
the problem in the coordinates with respect to the linear Hamiltonian Hq = —A + V, 

The nonlinear term ip'^tp, (assume A = 1) can be split into a sum of many terms using 
this decomposition. However, we claim that there is only one important nonlinear 
term in the equation for each component: 

ix = eox + (00, + • • • (1-18) 

= eii/ + (0i,2(x0o)(#i)O + --- (1-19) 
tdt^ = ifoe+ Pc''°a;|/Vo0? + --- (1.20) 

From ( |1.18| ), we know u(t) = e*^°*x(t) has less oscillation of lower order than x{t). 
Hence we say x{t) has a phase factor — cq. Similarly, y{t) has a phase factor — ei. The 
nonlinear term x?/^0o0i has a phase factor cq — 2ei, which, due to the assumption 
( p,.3D , is the only term in tp'^tp with a negative phase factor. It gives a term in ^: 

m = xy\t)<^ + ■ ■ ■ , $ = -— ^- 000? 

iJ + Co — 2ei — 

Notice that $ is complex. Substituting this term into ( |1.18D and ( |1.19| ), we have 

ix = i'jolyl^x + ■ ■ • , 
ly = -2«7o|x| \y\ y^ , 

with 7o given in (|1.8|) . In ( 1.21| ) we have omitted two types of irrelevant terms: 

1. Terms with same phase factors as x or y: for example, eox and \y\'^x in ( |1.18| ). 
Since their coefficients are real, they disappear when we consider the equations for 
|x| and \y\. 

2. Terms with different phase factors: for example, xy'^ in ( |1.18| ). Since these 
terms have different phases, their contribution averaging over time will be small. 
This can be made precise by the Poincare normal form. 
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From ( |1.21| ) we obtain the decay of y and the growth of x as well as the three 
time regimes mentioned previously. However, it should be warned that this set-up is 
only suitable when both x and y are of similar sizes. 



2 The initial layer and the transition regimes: The 
set up 

We now outline the basic strategy for the initial layer and the transition regimes. We 
first review the construction of the bound state families. 

2.1 Nonlinear bound states 

The basic properties of the ground state families can be summarized in the following 
lemma from 



Lemma 2.1 Suppose that —A + V satisfies the assumptions (AO) and (A2). Then 
there is an tiq sufficiently small such that for E between cq and cq + Xn^ there is a 
nonlinear ground states {Qe}e solving (|1.4|) . The nonlinear ground state Qe is real, 
local, smooth, X~^{E — Cq) > 0, and 

QE = n<j)o + 0{n') , n ^ C [X'^E - eo)]'^\ C = { j <PldxY^/\ 

Moreover, we have Re = OeQe = 0{n^'^)QE + 0{n) = 0{n^^) andd^^QE = 0{n~^). 
If we define ci = {Q,R)~^, then ci = 0(1) and Xci > 0. 

This lemma can be proved using standard perturbation argument and similar 
conclusions hold for excited states as well. For the purpose of this paper, we prefer 
to use the value m = Qi) as the parameter and refer to the family of excited 
states as Qi{m). It is straightforward to compute the leading corrections of Qi{m) 



via standard perturbation argument used in proving Lemma |2.1| . Thus we can write 
Qi as 

Qi{m) = m0i + (m^q^ + q^^\m)^ := m0i + q{m) , q^^\m) = 0{m^), q{m) ± 0i , 

(2.1) 

where gs = —X{Hq — ei)~^7r0f and vr is the projection 

Tih = h- {(t)i,h)<t)i. (2.2) 
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Similarly, we can also expand Ei{m) in m as 

Ei(m) = ei + ^1,2^2 + Ei,4m^ + E[^\m), E[^\m) = O(m^) . (2.3) 



Moreover, we can differentiate the relation of Qi{m) w.r.t. m to have 

'(m) _L 

(2.4) 



Q'i("^) = = 01 + 9'("^) , ?'(H = ^^("^) = 0{m^), q'{m) ± 0i, 



2.2 Equations 

Thus in the first and second time regimes, we write 

m = x(t)0o + Qi(m(t))e*®W + m , (2.5) 



where ^ E Hc(ifo)) see ( |1.12[ ). If we write 9(t) = 9{t) — Ei{m(s)) ds, we can write 
y{t) as 

y = me*® = mexp |i^(t) — i y _E'i(m(s)) c/s j> . (2.6) 

Denote the part orthogonal to 0i by /i = x0o + From the Schrodinger equation, 
we have h satisfies the equation 

idth = Hoh + G + A, 
G = X\^\^^ - XQle'^ 

= XQlie'^^h + 2h) + \Qi{e'^2hh + e-'^h'^) + X\h\'^h (2.7) 
A = [eq^ - ifnQ'^ e^® . (2.8) 



Since m(t) and 6{t) are chosen so that ( p.5| ) holds, we have = {(f)i,idth(t)) = 
(^01, G + {OQi - imQ[)e'^y Hence m(t) and ^(t) satisfy 

m= {(j)i,lmGe-'^) , ^ = --^ (0i, Re Ge"*®) . (2.9) 
We also have the equation for y: 

iy = irhe''® — {6 — _E'i(m))me*® = Ei{m)y + e*®(im — ni9) = Ei{ni)y + (0i, G) . 

Here we have used (^]9|). Denote A^ = ttA. We can decompose the equation for h 
into equations for x and Summarizing, we have the original Schrodinger equation 
is equivalent to 

ix = eoX+ (007 G + A^) 

ty = E,{m)y + {(P,,G) (2.10) 
tdt^ = Ho^+ Pc(G + AO 
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Clearly, x has an oscillation factor e and y has a factor e ^'^^^ Hence we 
define 

X = e-^^o*M, y = e-'^''v . 
Together with the integral form of the equation for ^, we have 

v = -te'^^' [{E^{m)-e^)y + {<P^,G)], 

^(t) = e-'^«*eo + / e-*^«(*-^) Pc^° Gi:{s) ds , G^ = t-\G + A,) 
Jo 

This is the system we shall study. 

2.3 Basic estimates and decompositions 

It is useful to decompose various terms according to orders in n so that we can identify 
their contributions. We now proceed to so this for G, A^, E{\y\) and ^{t). We expect 
that x,y = 0{n) and ^ = Oijn?) locally. 
1. G Recall that G is given by 

G = XQlie'^^h + 2h) + \Q^{e'^2hh + e^'^h'^) + X\h\'^h 

with h = x(j)Q + ^ and Qi = Qi(\y\). From the the decomposition (^TTI) of Qi = 

\y\(f)i + \yfq3 + Q^^Klyl), we decompose G as 

G = A(i/V? + 2y'y(l),qs) h + A(|i/| V? + 2|i/| Vi^a) 2/i 

+ A(#i + 7/2^93) 2|/i|' + A(#i + yy^gs) + A|/ip/i + (*) 

where (*) = A [2\y\4>iq^^^ + {\y\% + q^^^)^] {e'^^h + 2h) + Xq^^^ {e'^2hh + e~'®h^) with 
g(5) — g(^)(|y|). We then substitute h = x0o + obtain 

G = Gs + + G7, (2.15) 

where 

G3 = A(?/2x + 2|?/|2x)0o0? + A(2|xp?/ + x2y)0^0i + X\x\'^x(j)l (2.16) 

G5 = A(2t/3p + 4\y\^x) 0001^3 + A(2|x| + <Plq; (2.17) 
+ A(x0o + #i)'e + 2A|(a;0o + #i)|'e 
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(2.11) 

(2.12) 
(2.13) 

(2.14) 



and 

G, = A [2\y\M^'\\y\) + i\y\% + q^'\\y\)f] {e^'^'h + 2h) (2.18) 
+ Ag(^)(|y|)(e*®2/i^ + e-'®/i2) 

+ A (2y=^#ig3 C + 4|y| Viga + A(#i2|^r + #1^') 

+ Ay^ygg 2(a;0oe + ^^o^ + ICH + Ay^'^s (2x00^ + C') 

Note that G3 = 0{n^), = 0{rv') and G-j = 0{n^). If we use the convention 
that 

/ ^ 91 + 92^ 

for 11/11 < C \\gi\\ + \\g2\\ + ■ ■ • for some suitable norms, we have 

G < n^x + n^i + (2.19) 
G'5<n^x + n2^ (2.20) 
G7 < n^x + n^^ + n^^ + ^'^ (2.21) 

It is crucial to observe that no term in G3 is of order y^. This is due to our setup 
emphasizing the role of nonlinear excited states. The price we pay is the introduction 
of terms involving and g*^^^ . 

We now identify the main oscillation factors of various terms. For example, y^x = 
gi(-2ei+eo)t^2^^ and its factor is — 2ei + cq. For terms in G3 we have 

y'^x \y\^x \x\'^y x'^y \x\'^x ,^ ^2) 

-2ei + Co -Co -ei -2eo + ei -Cq 

From the spectral assumption |eo| > 2|ei|, — 2ei + eo is the only negative phase factor. 
Hence it is the only term of order (n^) that have resonance effect when we compute 
the main part of ^. Also, since has the same phase as it will be resonant in 
the y-equation. Similarly, \y\^x and \x\^x have same phase as x and will be resonant 
in x-equation. 

2. and E{m) Recall A^ = T:{eQi - imQ[)e'^. Since 9 = 0{n-^ \\G\loc) and 
^-0{\\G\U, 

\\K{s)\\ = 0(6) 0{7rQ,) + 0(m) 0{7rQ[) < Gn' \\G\\,^^ , (2.23) 
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To find out the main part of A^r, we substitute equation (|2.9|) for m and 6 to 
obtain ( m = \y\), 



= TT{eQi - imQ[) e^® 

= - {(01, G/2)m-'nQi + (01, G/2)m-i7rgie2^®} 
- t {(01, G/2t)7rQ[ + (01, G/2z)7rQ;e2*®} 



Since G = Gs + (Gs + Gy) and nQi{m) = m^q-i + q^^\m) by (PTTD, we liave ttQ^I 
3m^g3 + O(m^), and tlie main part of A,r is (also recall y = me*®) 

An,5 = -\ {(01, Ga) |i/|'g3 + (0i, G3) y'gs} 

-|{(0i, G3)3|i/pg3 + (0i, G3)3i/2g3} 
= -2g3(0i, G3|yp + G3y^) 



(2.24) 



Let Ajrj = A^ — A^^5. We have 

A, = A,,5 + A,,7 (2.25) 

< WGAoM? <n'x (2.26) 

< \\G, + G,l^M'+\\G\LM' ■ (2.27) 
3. ^ Recall the equation for ^ in (|2.14|) , 

^{t) = e-'^°*eo + /* e-'^°(*-^) Pc^° G^{s) ds , G5 = ^-^(G + A,) . 
Jo 

Since HAttH < Gn^ WGWi^^, the main terms in Gg = 'i^"'^(G + A^) is i^^G^. We now 
compute the first term Xy'^x(j)Q(j)l in G3 using integration by parts: 

-tX [ e-^^«(*-^) Pc?/^^0o0? ds 
Jo 

Jo 



—iXe 



-iHot ^ 



1 



i(i7o - Oi + eo - 2ei) 



y^x^i - e"*^°*2/^x(0)$ 



i(i^o -Oi + eo- 2ei) 
1 



giHos gj(eo-2ei)s 



(is 



Pc0o0icis 



g-iHo(i-s) gi(eo-2ei); 







(is 



(f ^m) $1 ds 



where 



$1 



-A 



Ho -Oi + eo- 2ei 



-cVOVl 



(2.28) 
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This term, with the phase factor cq — 2ei, is the only one in having a negative phase 
factor (see (|2.5|)). Since — (eg — 2ei) is in the continuous spectrum of Hq, HQ + eo — 2ei 
is not invertible, and needs a regularization —Oi. We choose —Oi, not +0i, so that 
the term e~*^''*?/^x(0)$i decays as t ^ oo. See Lemma p.2| . 

We can integrate all terms in and obtain the main terms of C,{t) as 



= y^X^i + |l/|^X$2 + \x\'^y^3 + X^?/$4 + |x|^X$5 (2.29) 



where 



— 2A — 2A 
'^'2 = -Tz Pc0o0i , $3 = T7 Pc0o0i , (2.30) 

-no — Co -no — Ci 

Ho - 2eo + ei Hq - eo 

The rest of ^ (t) is 

+ / e-'^»(*-^) Pe(G'5-z-^G'3-rU|e|'0 
io 

+ re-^o(*-) p^(riA|e|'0 
io 

= er(t) + d'\t) + (i) + ^4\t) + (t) 

where ef^(t) and ^P{t) are higher order terms in G^ which we did not integrate and 
the integrand G^ in ^^\t) consists of the remainders from the integration by parts: 

as ^ ds ^ 

as as as 

Here we single out C,^\t) since is a non-local term. Thus we have following 
decomposition for ^: 

m = ^^'\t) + e\t) = e(^) + (ef ) + . . . + ) , (2.32) 
2.4 Linear estimates 

We now summarize known results concerning the linear analysis. The decay estimate 
was contained in and |T^; the estimate ( |2.34|) was taken from and [I! 
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Lemma 2.2 (decay estimates for e **^'') For q E [2, oo] and q' = q/{q — 1), 



LI' 



(2.33) 



For smooth local functions and sufficiently large tq, we have 

1 



lim 

0-^0+ 



{Ho + Co - 2ei - ai 



<C{t) 



-9/8 



(2.34) 



L2 



where k = 1,2. 



Intuitively, we can write 

POO 

^^°+ Jo 

3 The initial layer and the transition regimes: The 
estimates 

In this section we wish to show the following picture for the solution ip{t): In the initial 
layer regime, the dispersive part gradually disperses away, while the sizes of the bound 
states do not change much. In the transition regime, the original dispersive part 
becomes negligible, while the 0o-components of ip{t) increases and the 0i-component 
decreases. 

Recall the orthogonal decomposition ip(t) = x0o+Z/0i+^ (|1.11|) . We have |x(t)p + 
bWP + IliWll^a = < n^- If we decompose ^ via ([LT2D 



I.e., 



^(t) =x0o + Qi(y)+e 



(3.1) 



we have y = y, x = x + 0{y^) and ^ = ^ + 0{y^). Thus 

|x(t)|, \y{t)l \m\\L-<ln, ||eo||y<4n 



For p = 1, 2, 4, define the space 



(3.2) 



-^loc 



{/: (x)-^^°/(a;)GL^(M3)} 



(3.3) 



where ro > 3 is the exponent appeared in the linear estimate (|2.34| ) in Lemma ^ 



Since ro > 3, we have 



-1 <C 



LP 



for any p, 



The following proposition is the main result for the dynamics in the initial layer 
and the transition regimes. 
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Proposition 3.1 Suppose that V satisfies the assumptions given in §i. Let %l){t,x) 
he a solution of ( |1.1| ) with the initial data ipo satisfying ( |1.14D . Let > be a 
sufficiently small constant to be fixed later. Let to = £3^"^. Then there exists ti and 
t2 such that for some constant C < 10000 we have 
1.01 , / n 



to < ii < 



Ian" 



log 



ti + C (nV)-i <t2<t 



10100 (7onV)^^ , 
(3.4) 



and the following estimates hold: 
(i) ForO<t< 2t2, 



\x{t)\ > I sup \x{s) 



0<s<t 



\m\\L^<C2nH'/'\x{t)\+C2Uo\\Y{tr'^' 
imilLf <C2n'\x{t)\ +C2Uo\\y{t)-'^' , 

loc 

\^^'\t)\\^, < C2n''/'\x{t)\+C2\My{t)-'^' ■ 

loc 



where the constant C2 will be specified in (|3.19|) of next subsection, 
(a) (Initial layer) for < t < to, 

l\xo\ < \x{t)\ < flxol, 
0.99\yo\<\y{t)\<1.01\yo\ 



1 /2 

(Hi) For ui = (|a;oP + | Iz/oP) defined in ( |1.16|) , 

> O.Oln , |x(t2)| > 0.99ni , \en < \y{t2)\ < 2en . 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Notice that (|3.4|) implies 



t2 < C^n 



(3.9) 



for some constant C3. 

We will prove these estimates using ( p. .141 ), ( |3.2|) and a continuity argument. Hence 
we can assume the following weaker estimates: For < t < 2t2- 

\x{t) \ > I sup |a;(s)| , 

0<s<t 

|x(t)| < 2|xo| for t <to , 
\\m\\L^<2C2nh'/'\xit)\ + 2C2\My{t)-'^' , 
\m\\Lf <^C2n'\x{t)\ +2C2Uo\\Y{tr'^' , 

loc 

U^'\t)\\^, <2C2n''/'\xit)\+2C2\MY{tr'^' ■ 
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By continuity, if we prove Proposition |3.1j assuming these weaker estimates, we have 



proved the proposition itself. We shall see also estimates will be used only in 
estimating higher order terms. 

Recall from (|2.23| ) that the local term satisfies HAttH,^ < Cn"^ \\G\\^i for any r. 



Thus we have 

m\ < \\Ghi > m\ < \\Gh. +\y\'. (3.10) 

loc loc 

The following lemma provides estimates for G assuming the estimate (01). 



Lemma 3.2 Let G be given by ( |2.15| ) -( p. 181 ). Suppose n is sufficiently small and the 



estimate (^ holds for t < G^n ^. Then we have the following estimates for G: 

||G(t)||^4/3nL8/r < G,n'\x{t)\ + GiG^y UolW (ty'^' ■ (3.11) 

ll(G'-G3)(t)|L4/3nL8/. <C4ni^/1x(t)|+C(C2)n2||eo||y(t)-'/' • (3.12) 
\\{G - G.mWL} <C,n'\xit)\+G{G,y\MY{tr'^' ■ (3.13) 

loc 

where C4 is a constant independent of G2 and G{G2) denotes constants depending on 

C2. 

Moreover, (|3.11|) and ( |3.12|) remain true if we replace G by G^, and {G — G3) by 



{G^ — i ^Gs). Furthermore, we have 

\\G^mL^<C,n' ■ (3-14) 

By the assumption ( |1.14| ), when t > to the last term Gn^ ll^olly {t)~^^^ is smaller 
and can be removed. The proof of this lemma is a straightforward application of the 
Holder and Schwarz inequalities. 

We will use (^) and (PTT^ ) for ^, and ( [XT^ ) for x and y. 
Proof. Recall ( |2.15| ) that G = G3 + G^ + Gj. We first consider only the non-local 
term A|^p.^ in G, which belongs to G-j. Since n < and ^2 ^ C'^e^'^n^'^, by (^) 
we have ||e(s)||^4 < Gn^/^\x{s) \ + G ||^o||y {s)~^'^- Also using (jSj) and the Holder 
inequality we have 

\\m{s)\l^^. < C \m\\% < C {n'/'\x{s)\f + G Uofy {s)-'^' • (3.15) 
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Li 

<Cn'l^ {[n^l'\x{s)\f^ W'''') (3.16) 



Hence this non-local term satisfies (|3.11|) -( ^.12D . Moreover, to prove ( p.l3| ), we can 
bound|||ePe(^)||.i by 11^(^)11^4. 

loc 

For the local term G— A|^p^ = G3+G'5 + (G7 — A|^p^), all L^-norms are equivalent. 
We can read from the explicit expressions of G the following estimates: 

G, < n'^x + nH 
Gi-m'i < n'\x\+nH + ne 

To estimate ^ in — A|,^p^, we can use W^W^i . For example, 

loc 

ll^ie'lLva < C\y\ m% < Cn ({uM? + (s)-'^' 



Together with the explicit expressions of G and G3, similar arguments show ( 3.11 )- 



Since G^ = i-^iG + A^) and ||A^|| < Gn^ WG^^^ by QZ^ ), (|XTT| ) and ( PT^j) hold 



if we replace G and G — G3 by G^ and G^ — i ^G^. obtain ( p. 14] ), we only need to 
check the non-local term A|^p,^. Since ||^(s)||^4 < {2G2GI^^ + 8^2)^, we have 



lliereiLi < C4,i ms)h. imwl < ^c,nn' , (s.i?) 

provided we choose G5 > 10(2^2^3''"' + 8^2)^4^!. Thus the lemma is proved. Q.E.D. 



3.1 Estimates of the dispersive part 



We now prove the estimates for in Proposition |3.1| by using (p.2|), (H), Lemma 2.2 



and Lemma 3.2 



Step 1. and Lf^^ norms, < t < 2t2 
Recall the equation ( |2.14| ) for ^: 



^(t) = e-'"°'io + / e-*^«(*-^) Pc^° G^{s) ds , G^ = r^G + A,) . (3.18) 
^0 
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By (PD , Lemma and Lemma ^.2|, we have 



Ile(t)|lz.4<||e-^^%|L,+ / C\t-sr'''\\G^{s)\Uds 



-3/4 



3/4 



4/3 



1-3/4 



C,n'\x\is) + CiC2)n'Uo\\Y{s) 



'9/8 



ds 



<C,,, \My{t)-"' + C2,inH'l'\x{t)\ + C{C,)n'\My{tr'''^ 



where 6*2,1 is some exphcit constant. 

We now estimate H'Cl^)!!/,* . If t < 1, we can bound Lf^^-norm by L^-norm. Hence 

loc 

we may assume t > 1. For , We divide the time integral in (|3.18|) into s G [0,/!: — 1] 
and s E [t — 1, t]. In the first interval, we estimate Lf^^ by norm; in the second we 
estimate Lf^^ simply by L^. Using similar arguments in the previous estimate of the 
L4 norm, we have 



loc 

<C\My{t)-"' + 



t-1 



U _ _5|9/8 ll'^dlLS/y ds + 



_^ \t - s|3/4 ll'^dlLVS ds 



C 



{t-s) 



9/8 



C^nM^s) + C^C^W UoWy {sy'"') ds 



+ sup 1 1^5(5) 11^4/3 

t-l<s<t 



where C2,2 is some explicit constant. 
Step 2. Lf„^-norm, < t < 2t2 



Recall the decomposition (pl): ^ = ^(2) + ^0) = ^(2) + ^^(3) + . . . + ^(^)). We will 



estimate the L^^^-norm of each term. 



0. Since $1 G L^^^, and <l>j G L^, (j > 1), we have 

U^'\t)h. < C^^sCn' \x{t)\ , 

loc 

for some explicit constant 6*2,3. 

1. We have 

loc 



for some explicit constant 6*2,4 by the L^''^ estimate of e in Lemma p. 

2. ,^2'^''. By the linear estimate (p.34|) in Lemma |2.2| we have, for some constant 

^2,5, 
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-9/8 



loc 



3. To estimate ^i^\t) = - /J e-'^"^^''^ PcG^ ds with d defined in (|23T|) , we 
need estimates ( |3.1U| ) for ii and v and tlie linear estimate ( p^.34D in Lemma pj.2| . Hence 



< / ||e-'^«(*-")PcG4||^2 ds 

1 loc 



LL Jo 

3^ C / (t — s)^^^^ (n^l-itl + n|x?)|) (is 

(3 10) 

V C / (t - s)-^/^ {n^ ||G'||^4/3 + nVl) ds 
Jo 

& ^ f\. .-9/8/ 4, , , ^.^^ 4„. „ 



^c}\t-sr^^{. 
<C2,en''\x{t)\ + C{C2)n'\MY{t) 



-9/8 



4. d^^ + ^f^. We write + ^^"^ = Joe-'""^''' 

i^^Os){s). By Lemma p.2| and Lemma |2]^, we liave for t > 1 



ds, where G^^^^s) :- 



< 



.0 .v^. 

loc 



(is 



<o O 



\t - s\ 11^5,5(5)118/7 ds + 

Jt-i 



<o(r\t-s\-"' 

\J0 Ul~L 

<02,7n''''\x{t)\+0{02WmY{t) 



Jt-l 

„^^2 ||. |, _ /,.-9/8 



'*^C|t-s|-^/^||G5,5(s)||,/3 ds 



'''''\x{s)\ + 0{02)n^\M (5)"'/') ds 



for some exphcit constant 6*2,7. If t < 1, we can bound the Lf^^-norm by the L^-norm. 
Hence the last estimate for t <1 follows from the estimate in Step 1. 

We have obtained estimates on ^ involving explicit constants 6*2,1, • • • ,6*2,7 and 
C*(C*2). We now define the constant O2 in (^]6|) to be: 



O2 = C*2 1 + ■ ■ ■ + C*; 



2,7 • 



(3.19) 



Since all terms involving C*(C*2) have some extra n factor, ^ satisfies the estimates in 
( ^.61 ) provided that n is sufficiently small. 

Summarizing, we have proved the following lemma: 

Lemma 3.3 // n is sufficiently small, there is an explicit constant O2 such that, if 
(I) holds in [0,t] for some t < O^n^'^, then the estimates (|3.6| ) in Proposition \3. 1\ also 
hold in [0, t]. 
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Remark (3.20) In the proof, we only used ( |3.2| ), Lemma and Lemma 
73i. The information we need on the size of bound states is in (13.21) and the first 



estimate of (|). Since (|3.2| ) is always true, we only need to ensure that the first 
estimate in (||) holds. 

3.2 Normal form for equations of bound states 

We now compute the Poincare normal form for the bound states. This normal form 
will be used to estimate the bound states components x and y in next subsection. 
Recall that we write 

x{t) = e-'^^'uit), y{t) = e-''''v{t) 

and the equations ( p.l2|) and (|2.13| ) for u and v, 



u = -le''^' (00, G + A^), v = -ze'^i* [(Ei(m) - ei)y + G)] 



Using the decompositions ( p.25| ) for A^r and ( |2.3| ) for Ei{m), we can decompose the 
equations for u and v according to orders in n: 

= -R«,3 + -Rm,5 + RuJ , (3-21) 
-^e-^*[(0i,G,) + Ef)(|i/|)y] 

= Rv,3 + Rv,5 + Rv,7 ■ (3.22) 



Using ( p.26| ), ( p.27| ) and Lemma [3^, which assume (H), we have 



Rujl 


< 


WGr 


Rv,5\ 


< 


\\G5 


Rvjl 


< 


WGr 



\Ru,5\ < 11^511^1 +||A.,5L, < n'\x\+n'my{t)-'/' (3.23) 

loc 

,^ + II A.jLi < n^\x\ + UoWy (ty'^' (3-24) 

+ \y\' < + UoWy {ty'^' (3.25) 

'C 

< WGAlI + \y\' < + \My it)-'/' (3.26) 

loc 

We shall first integrate Ru,3 and R^^s in step 1, and then integrate Ru,5 and R^^^ 
in step 2. 

Step 1 Integration of terms of order 

In the equation of u, ( p.21| ), the terms of order are contained in i?„3 = 
— 'ie*^°*((/)o, Gs). The resonant terms from G3 are \y\'^x and whose phase cancels 
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the factor e*^°*. The other four terms of order in have different frequencies and 
can be exploited using integration by parts. By ( p^.l6D we have 



where 



= ci + C2 \v\'^u + ^(^r) + 9u,l + Ru,5 + Ru; 



3i{2eo-2ei)t 2-, j(eo-ei)io| 12 j(-eo+ei)t 2 



= - A0O, 7^ 0001 H 0001 H ^ 

\ 2eo - 2ei cq - ei -cq + ei 

and 



i(2eo-2ei)t d/ 2^N e^^""""!)* ^ (2 |m Pt;) „ e*(-"«+"i)* ^ (w^t;) 
= A0O, ^ ^000? + ' ' V o01 + -^0^1 



2eo - 2ei cq - ei -cq + ei 



In the equation of v, ( p.22|) , the terms of order are in i?^, 3 = — 'ie*^i*[(0i, G3) + 



-E'i,2|l/Pl/]- There is only one resonant term in Gs, namely, \x\'^y. Another resonant 
term of order ra^ is from the term Ei^2\y\'^y- The other four terms of order in G3 
have different frequencies and can be integrated. We thus have 

V = -2e^^^*[(0i, G3) + E,,2\y\%] + Rv,5 + Rv,7 

= ce\u\'^v + cr\v\'^v + ^(^f) + 9v,i + Rv,5 + Rvj 

where 

C6 = -i2A(0o, 0i), C7 = -iEi^2 



^i{-ei+eo)t„.2^ „i{ei 



= - \ A01, ■ 0001 + 



V u , ,0 e^^ z\v\ u 



eo)t2|^ 



ei + eo ei - cq 



gi(2ei-2eo)t^2^ gi(ei -eo)t ^ 1 2^ 

+ ^ 0001 + ^0 

2ei — 2eo ei — eo 



and 



= A01, ^000? + ' ' ' 000^ 



-ei + eo ei - eo 

+ ^ f^0 01 + ^^^^0 

2ei — 2eo ei — eo 
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We now define 

Ui = u — Ui , Vi = V — . (3.27) 
The equations for ui and Vi are 

ill = Ci\u\'^U + C2\v\'^U + Qu,! + Ru,b + Ru,7 

= C]\ui\^Ui + C2\V]\'^Ui + gu,2 + gu,l + -R«,5 + Ru,7 

gu,2 = Ci(|m|\ - + C2(|t^|^M - jt'll^Ml) 

and 

Vi = C6\u\'^V + C7\v\'^V + Qy^i + i?^,5 + Ryj 

= Cel^il^Vl + C-j\vi\^Vi + fl'^^a + + Rv,b + -Ri;,7 

5'i;,2 = Cq{\u\^V - |ltl| V) + C-j{\v\^V - |fl| V) ■ 

We have finished the integration of order terms. Note that both u'^ and v'^ 
enter the equations of Ui and v^. This is the reason we compute their normal form 
together. 

Observe that 



^ nM , \v^\ < n>| . (3.28) 



We now decompose g^^i, Qv,!, gu,2 and g^^2 according to their orders in n. We want to 
write them as sum of order and order terms. We first claim that gu,i and g^^i 
are of the forms 



9u,l — e**^°*fl'u,l,5 + 9u,l,7 1 9v,l — ^^^^'^9v,l,5 + 9v,l,7 , 

where gu,i,7 and 5^^,1,7 are order -nJ terms, and gu,i,5 and gv,i,5 are explicit homogeneous 
polynomials of degree 5 in x, x, y, y with purely imaginary coefficients. Moreover, 
every term in gu,i,5 has a factor x or x. For example, the first term in gu,i is 

= C'e^(2-o-2-i)*(2wvv + v'^u) 

= Ce''°' (2xye-'^i*[i?,,3 + Rv,5 + Rv,7] + y'e^'^^l^ + ^ + ^]) 
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where C = (A0o, (2eo — 2ei)~^0o0i) is real. Repeating this calculations for all terms 
in (yf„ 1 and collecting terms of order n^, we obtain gu,i,5- The rest belongs to gu,i,7- 
There are two terms of order in the last expression: the terms of the form 2xyRy 3 



and y'^Ru^s. By definitions of i?„,3 and Rv,3, they are explicit polynomials of degree 5 
in X, X, y, y with purely imaginary coefficients. 

From the estimates of ( p.23|) - (|3.26| ), we can bound gu,ij by 

k,i.7(t)| < n'\u\ + n'\My{tr'^' , ^ ^ 

(o.29j 

Similarly, we can write gu,2 and (7^,,2 as 

gu,2 = e^^°^gu,2,5 + gu,2,7 , gv,2 = ^^^^^gv,2,b + gv,2J , 

where gu,i,5 and gv,i,5 are explicit homogeneous polynomials of degree 5 in x,x,y,y 
with purely imaginary coefficients and gu,2,7 and gv,2,7 are order rJ terms satisfying 

\9u,2,7{t)\ < n^\u\ , \g.,2,7{t)\ ^ Av\ ■ (3.30) 

Here we have used (|3.28|) in last estimate. Moreover, every term in gu,2,5 has a factor 
X or X. We shall not perform calculations and estimates in details as they are similar 
to the previous step. To gain some idea, we shall do one example and show it is of 
the right form. By using u — Ui = , the first term in gu^2 can be written as 



u\ u— \ui\ Ui = u u — {u — Ui) {u — Ui) = u (mi ) + 2\u\ Ui + 0{u\ui I ) 



x^{u^ ) + 2\x\'^u^ + 0{u\u^ p) 



The first two terms, x^(m^) + 2|xpMj~, contributes to 5^^,2,5 • Since equals to e^^°^ 
times a polynomial of degree 3 in x, x, y, y with real coefficients and Ci in 2 is purely 
imaginary, they are of the desired form. 
Summarizing, we can write 

gu,i + gu,2 = e*^°*-Ru,5 + gu,3 



gv,i + gv,2 = e*'i*i?„,5 + g 



v,3 



where Ru,5 = 5 +fi'«,2,5 and Ry^^ = gu,i,5+gu,2,5 are explicit homogeneous polynomi- 
als of degree 5 in x, x, y and y with purely imaginary coefficients. Moreover, every 
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term in Ru,5 has a factor x or x. Also, gu,3 = gu,i,7 + gu,2j and 5(^,3 = 5(^,1,5 + g^^2,b- 
From the assumption (|^), we have 

\RuA < ^Vl, |^,,5| < n\ (3.31) 

|^?n,3| < n^\x\ + lleolly {t)-'/' , \g,,s\ <n' + n' \My {t)-'^' . (3.32) 
The final equations for ui and f 1 are 

ill = Ci\ui\\i + C2\vi\\i + {Ru,5 + e'^°'Ru,5) + {Ruj + gu,3) (3.33) 
7)1 = ce\ui\%i + crlvil'^vi + (i?„,5 + e'^^'R^^^) + {R^j + c/^,3) (3.34) 

Step 2 Integration of terms of order 

We now integrate terms of order n^. In wi-equation (|3.33|) we have Ru,5 + e^'^°^Ru,5, 
where Ru^^ is from the decomposition of original equation ( p.21|) and Ru,5 is from the 



error terms gu,i + gu,2- Similarly, terms of order in fi-equation (|3.34|) is /2„^5 + 
e^'^^^Rru^^. Observe that they are either of the form x'^y^ with |a| + \f3\ = 5, or of the 
form xyC,. Also note that there are two sources in Ru,5- G5 and A^^s. Among all these 
terms the main term is G5. 

We have already studied -R„_5 and Rv,5- They are explicit homogeneous polynomi- 
als of degree 5 in x, x, y and y with purely imaginary coefficients. Moreover, every 
term in Ru,5 has a factor x or x. 

We next look at A^. Recall (|2.25| ) that A^ = A^^s+A^j and A^^s = -2^3 (0i, G3\y\'^+ 
G^y"^) ( |2.24| ). Thus A^r^s is a homogeneous polynomial in x, x, y and y of degree 5 with 



purely real functions as coefficients. Therefore the term — ie*^°*((/)o, A^.5) in ii equation 
( ^.21| ) gives only polynomials with purely imaginary coefficients and a phase e*^°*. 
Recall G5 is given by 

G5 = X{2y^yx + 4||/|^x) 0001^3 + K'^lxWy + x'yf) (gq^ 
+ A(x0o + i/0i)'e + 2A|(x0o + l/0i)|'e 

Recall the decomposition ^ = + where 

i^'^\t) = y'^x^i + \y\'^x(^2 + \x\'^y^3 + x'^y(^4 + |xpx$5 

with $1 the only function with nonzero imaginary part ( p.28| ), ( p.3(]| ). Let 

$1 = $1 + i^^j (3.35) 
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with both $1^^ and real. Denote the real part of ^^'^\t) by 

We can write ^ = y'^xi^j + ^^^^ + Thus we can further decompose G5 as 



(3.36) 



G5 — Gs,! + (^5,2 + G^^s 



(3.37) 



where 



^5,1 = (x0o + y(l)i)Yxi-i)^i,i + 2|(x0o + #i)| V^^'^'i,/ 
6-5,2 = A(2?/3p + 4|y|^a;) ^o^igg + A(2|x| + a^'z/y") 0o?3 

+ A(x0o + l/0i)'d? + 2A|(x0o + #i)|'el?^ 
G5,3 = A(x0o + #1)'^ + 2A|(x0o + #i)|'e^'^ 

The term 3 will be shown to be smaller than 6*5,1 and 65,2- Although 65,1 
and 65,2 are of the same size, G5.2 consists of monomials in x, x, y and y with real 
functions as coefficients, while 65,1 with purely imaginary coefficients. The reason 
that 65,1 has purely imaginary coefficients is due to the resonance of some linear 
combination of eigenvalues with the continuum spectrum of Hq appearing in the 
form {Hq -Oi- 2ei + eo)~^ 

The only resonant term in w-equation from 65,1 is \y\'^x (from y'^C,)'- 

-ze*^o*(0o, (#1) Va:(-0$i,.) = -{(Po<Pl^i,i)\v\'u , 
and the only resonant term in f-equation from 65,1 is (from xy^): 

-ie'''\^l,2ix(|)o)iy(pl)y^x^^l,^) = 2{M^^l,^)\u\M'v • 



Note their coefficients only differ by a factor —2. We recall 

-A 



7o 



-(0001, $m) = -Im A0O0I, 



c<PO<Pl 



Hq — Oi — 2ei + Co 
Together with the definitions of Ru,5 and i?„,5 in ( |3.21| ), we can rewrite 



> 



(3.38) 



Ru,5 — i{4>0, ^5,1 + ^5,2 + A^,i 



2e^^°*(0o, ^5,3) 



i?„,5 + e*'^*i?.,5 = e'''' Rv,5 - ^(01, 65,1 + 65,2 + E,,4\yry) 



As in Step 1, we now integrate by parts the non-resonant terms inside the square 
brackets. The resonant terms can't be integrated and we shall only collect them. 
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This procedure is the same as in Step 1 and we only summarize the conclusion: there 
exists constants 03,04,05,03,09,010, U2 — 0{u^ + uv"^) and = 0{u^ + uv"^) two 
homogeneous polynomials in u and v of degree 5, and g^^^ and g^^A the integration 
remainders such that 

^,;,5 + e'^^*^„,5 = (csIhI^ + Ogj-Ul^l^ + Oio|v|^) V 

Furthermore, except 03 and 09, all other constants are purely imaginary. The real 
parts of 03 and 09 are from G^^i and they are given explicitly by 

R-eo3 = 7o, Reog = -270 . (3.39) 

The explicit forms of U2 or V2 are not important. We only need to know their sizes. 
We can now write the equations for u and v as 

■ui = oi + 02 + (o3|ii|^ + 04|ii|^|i'|^ + Osl^;!^) u 

+ ^ ("2 ) + 9uA - ie'^"\(po, ^5,3) + gu,3 + Ru,7 , 
Vi = 06 \uifvi + O7 \vifvi + {cs\u\'^ + Ogl^l^lvl^ + Oiol^l^) V 

+ ^ (^^2") + 9v,4 - ie^^i*(0i, ^5,3) + ^^,3 + R^j , 

We now define 

IJ, = ui — U2 — u — Ui — U2 (3.40) 

u = vi — V2 = V — Vi — V2 (3-41) 

We have 

/i = Ci l^il^Mi + C2 l^ll^Ml + (c3|m|^ + C4|M|^|t'|^ + Csl^l^) U 
+ Qua - ie'^°\(t>0, 05,3) + gu,3 + Ru,7 
= Oi |//| V + C2 V + (csIa*!^ + 04|/ini^|^ + Osli^l^) fJ' + gu 



and 



= Cg \Ui\'^Vi + C7 \Vi\'^Vi + (csImI^ + Cgl-Upl^P + Cio|t'|'^) V 

+ gvA - ie'^^^h, (^5,3) + gv,3 + Rv,7 
= 06 l/il'^iy + 07 + {cs\ii\^ + 09I//I + cio|i^|^) J^ + Qv 
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with 



Qu = gu,4 + gu,5 + gu,3 + Ruj - ie'^° (0o, ^5,3 



5,3) 



(3.42) 
(3.43) 



and 



9u,5 = Cl (IWlpMl - |/U|V) + C2 {\Vi\'^Ui - |Z/|V) 

+ Csd^l^M - + Qd^nt^pM - |/ip|z^P/i) + Csdf I^U - |z^|^/i) 

fi'^.s = C6 (Iwil^t^i - |/i|V) + C7 (|t;ipi;i - |z^|V) 

+ csdul^f - |/i|'^z^) + cgduplf l^t; - |/i|^|z/pz/) + Cio(|f l^'t^ - l^^l^z^) • 

Observe that the error terms Qu and Qy are of the form 

gu, gv ~ (a;' + Xy^) + (x^ + y^) ^(3) + + ^4)^ + ^3) ^ . . . ^ 

where denotes some local function. For we should add \y\^: Since ^f^, has a term 
_^g«eit ^0^^ £'(^)(|y|)y) from -Ri,j. This term has no factor in x and is of order \y\'^ . 

Under the assumption of (|^), the error terms gu^ and (^^^4 can be estimated simi- 
larly as for Qu^i and g^^i. Also, gu,^ and (71,^5 can be estimated similarly as for gu,2 and 
gv,2- We also have, for j = 1, 2, 



Ae^t I 



G,,s)\<Cn'\\e^\ 



^loc 



Cn6-^/>|+Cn2||eo||y(t) 



-9/8 



Together with the estimates (p.24| ), (p.26| ) and ( p.32D , we conclude 



\gu{t)\ < C,n'-'/'\xit)\ + C,n' UoWy (t)"'/' 
\g.{t)\ < Cen'-'/' + Cen' UoWy {t)-'^' 

Summarizing our effort, we have obtained the following lemma. 



Lemma 3.4 Let ji and v he defined as in ( p.40| ) -( |3.41| ). They satisfy 
fi = (cilyup + C2|z/H/i + (csl/il"^ + C4|/i| VP + cskH/i + gu 

1> = (Cel/ip + C7|z/p)z/ + (Csl/il"^ + Cgl/i] VP + CioW\^)l' + gv 

All coefficients Ci, ■ ■ ■ , Cio except C3 and Cg are purely imaginary and we have (p.39| ), 



(3.44) 



I.e., 



Re C3 = 7o 



Re Cg = -270 



(3.45) 
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where 70 > zs defined in 

Moreover, assuming (^) and using the estimates ( ^.21 ) and Lemma |J. 4 we /iai'e 



\guit)\ < C,n'^'/'\x{t)\ + C,n' UoWy (t)-'^' 
\g.{t)\ < C,n'~'/' + c,n' \My {ty'^' 



(3.46) 
(3.47) 
(3.48) 



for some explicit constant Cq. 

3.3 Estimates for bound states 

In this subsection we will conclude the estimates for x and y stated in Proposition 
3A\ . Recall tha, t under the assumption of (^ and (|3.2|) , we have proved Lemmas 
3l^- p.4| which contains estimates for ^, Qu and g^. 

We now derive some preliminary estimates. Let / = 2|/ip and g = |z/p. We have, 
by (pD, 



\f{t) - 2\x\'\ < 5Cen^\x\', \g{t) - \y\'\ < 5C,n' 



(3.49) 



We also have from ( p.44|) that 



/ = Reififi = 2jog^f + Reifigu , 
g = Re 2i/z> = —2'jofg'^ + Re 2z/(yf^ 



(3.50) 
(3.51) 



By ( 1^1 ) and (g), we have 



f + 9 



|Re4/2^„ + ReAug,\ < 4^6^^-^/^ + AC^^n^ UoWy (^) 



-9/8 



Recall nl = \x{0)\^ + l\y{0)\^. By ( |3^ we have |(/ + ^)(0) - 2nl\ < lOCgn^ Thus, 
for t < C3n-\ 

\{f + g){t)~2nl\ < WC,n''+ [Ucn'^-'/^ + ACen^UoWyisr^/Us 

Jo 

< ^CsCen^-^/^ . (3.52) 



We now prove Proposition 3.1 in three steps. 



1. Initial layer regime 

In this period the dispersive part disperses away so much that it becomes negligible 
locally. The times it takes for this to happen is of order = ^3^~^- We first prove 
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that ||xo| < \x{t)\ < ^\xo\ for t G [0,to]- The main ingredients of the proof are the 
norm form equation of x from the last section and the foHowing observation. The ^ 
dependent term is of the form n^^ or of higher orders. Because of our assumption 
IKolly ^ £2\xo\n''^ and the decay of ^(t), this term will not change x{t) very much. 
More precisely, for t G [0,to], = £3^"^, we have by ( g30D , ( gTzD , (|]4|), the 
assumptions ||^o||y < £2\xo\n~^ and (|^), 

|/(t)-/(0)| < r47o//(s) + 5|a;o||^7„(s)|rfs 
Jo 

< (87on^ + Cn"^-^/^) Ixol'to + lO\xo\n^ UoWy 

< (IO70 + 1)^2 |xo|' + lOeskoP < |/(0), 

provided n, and £3 are sufficiently small. By ( p.49|) , we have ||a;(t)p — \xq^\ < 
^|xop. Hence we have ^\xo\ < \x(t)\ < ||xo| for t G [0,to]- 

Similarly, we can show \g(t) — g{0)\ < ((IO70 + l)e2 + 10e3)n^. Hence we have 
\\y(t)\ — \yo\\ < 0.01|?/o| for t G [0, to] if £2 and £3 are small. The smallness of 62 and 
£3 can be guaranteed if we define 

_ 1 _ 1 

~ 2000(70+1)' ~ 2000 ■ ^ ' ^ 

2. Transition regime (i) 

In this period most mass of the disperse wave is far away and has no effect on the 
local dynamics; the ground state begins to grow exponentially until it has the order 
ra/lOO. The time it takes is of order n~^. 

Define 



We want to show that 



h = inf {t : \x{t)\ > O.Oln} := t[ . (3.54) 

t>to 



0<ti<to + ^ fogf^) , (3.55) 



Suppose ( p35D fails, that is, \x{t)\ < O.Oln for all t < t[. By (|^ and j ^A^ ), we 
have f{t) < 0.0004 and g{t) > 0.9995 ra^ for t < t[. Hence 

fit) > 270(0.9995 n2)V + 0(^6)/ > Y^7onV , 

if n is sufficiently small. Hence 

/(t) >/(0)exp|^7on^t| , (3.56) 
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for t <t[. We have 

fit[) > /(0)exp|-^7on^ ^ log ) 1^ = f{OW\xo\-' > 0.99n^ 

which is a contradiction to the assumption that < O.Oln. This shows that 

ti satisfies ( ^.551 ). We also have that ( |3.56[ ) holds for all t < ti, and that /(ti) > 
5 ■ 10-5 n^. 

3. Transition regime (ii) 

Recall the definition of e ( |1.15|) . Define 

t2 = inf {t : g{t) < {enf} . (3.57) 

We want to show that 

ti<t2<ti + 10100 (7oriV)"^ := 4 • (3.58) 

Suppose the contrary, then g{t) > (eri)^ for all t < tg. Then / > and we have 
that /(t) > 5 ■ 10-5 for ti < t < t^. Hence 

g < -(1.99)70/^?' < -9.95 ■ 10-^70^2/ • 

Hence 

git) < [gih)-' + 9.95 ■ 10-'^on\t - t,)]-' , (h < t < Q. 

and (7(^2) < (sny. This contradiction shows the existence of ^2 satisfying ( p.58| ). 

Since g > — {2.01) fg"^, similar argument shows ^2 > Ce^'^n"'^ if \y{0)\ > 2en. 
Combining with estimate ( |3.52| ) for f + g, we have estimates for /(t2)- From ( |3.49| ), 
these estimates of / and g can be translated into estimates of x{t2) and ^(^2) stated 
in Proposition ^.1| . 

We have proved (|3.4|) , ( ^^ , (|3^ ) and ( p. 81) in Proposition ^J] , using the estimates 
( p..l4| ), ( p. 21) and the assumption (|]). Since (|^) holds for t = 0, by continuity it holds 
for all t < t2. From Lemmas |3.2| - |3.4| and the above estimates, Proposition |3.1| is 
proved. 

4 Stabilization regime 

In this section we study the solution il){t) in the third time regime, after the solution 
has become near some nonlinear ground state. In this regime, it is natural to use the 
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decomposition ( [4 .31 ) for the solution ■ip(t) which emphasizes nonhnear ground states. 
A key issue here is to pass the information from the coordinates system ( |1.12| ) to the 
( [4.31) . As emphasized in the introduction, it is not sufficient to use only the estimates 
of ip at t = t2. We will also use the explicit form of the main terms in the dispersive 
part of tlj{t) to ensure that they do not come back to affect the local dynamics, i.e., the 
part of the wave represented by these terms is "out-going'. The set-up and proof here 
are similar to those in except the big terms of the dispersive part just mentioned. 



We shall first show that certain local estimates used in [|12| are still small (see Lemma 



4.1 Preliminaries 



In the time regime t > t2, we will use the set-up in [0, which we now briefiy recall. 
For the proofs we refer the reader to |12 . 



For all nonlinear ground state Qe with frequency E, Let £e be the linearized 
operator around Qe- 

Ch = -i{{-A + V - E + 2XQl)h + XQlh} (4.1) 

With respect to Ce, we can decompose L^(]R^, C), as a real vector space, as the direct 
sum of three invariant subspaces: 

L\R^C) = S{Ce)(BEi{Ce) © M^e) (4.2) 

where S{Ce) and Ei{Ce) are generalized eigenspaces, obtained from perturbation of 
00 and (pi respectively and Hc(£_b) corresponds to the continuous spectrum of Ce- 
Notice that this decomposition is not orthogonal. Also, S{Ce) = span^{iQE, Re), 
where Re = OeQe- 

For each ip sufficiently close to Qe, we can decompose ip as 

4^ = [Qe + aERE + Ce + Ve] e'®^ . (4.3) 

Here aE,QE £ and (e ^ Ei(£) and rjE G Hc(£). The direction iQE is implicitly 
given in Qe{g^^ — !)• Moreover, for this ip there is a unique frequency E' such that in 
the decomposition ( [4.3|) the coefficient a vanishes. In some sense it means that Qe' 
is the closest nonlinear ground states to ip. 
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4.2 Estimates 



Our aim is to show that ip{t2) satisfies the conditions of Theorem 3 (resonance dom- 
inated solutions) in ||12|. By Proposition |3.1| one can show that 'ip{t2) is close to a non- 



linear ground state QeoC*®" inL^^^-norm, i.e., HQeoIIls = '"'t2 < ""-o, ll'^ih) - e*®" 1 1 r 2 < 

loc 

£0^*2- From the same Proposition, we have rit^ ~ ni ~ n where ni is defined in ( |1.16|) 



and n is defined in ( |1.14| ). Thus for the purpose of order of magnitude, we are free 
to interchange nt^ with n. We now state the conditions of Theorem 3 (resonance 
dominated solutions) in [Q: Suppose that the initial data 4'(t2) is decomposed as in 
( [4.3| ) with the frequency E = Eq chosen so that the coefficient a vanishes. Then the 
excited state component ( satisfies 

< IICII < eon (4.4) 

and the dispersive part satisfies 

\\VE{t2)\\y < C WCf (4.5) 

for all E close to Eq with \E - Eo\ < \\C\\ . Here the Y norm is defined in ( [1.7|) . We 
shall see that the condition ( [4.4|) is easy to verify by Proposition The dispersive 
part, however, is no longer localized and there is no hope to satisfy ( [4.5|) . In fact, 
even its L^-norm is not small enough. Recall we know from the previous sections 
that |a;(t2)P — l^ol"^ is roughly one half of |?/oP ~ |l/(^2)P- Thus, by conservation of 
L^-norm, the dispersive part gains the other half of |?/oP ~ |l/('^2)P and 11^(^2)11 ~ n. 

One believes on physics ground that most mass of the dispersive part is far away 
and it has little infiuence to the local dynamics. The local part of the dispersive 
component, on the other hand, is generated by changes of the bound states and is 
small. Thus the results in should still hold. This idea, however, requires to 
clarify the concept of "out-going waves" for nonlinear equations. Instead of directly 



approach this problem, we examine the condition ([4.5|) in the proof in ||T2| . It is used 
only to guarantee that for all s > 

II sC /J. Nil ^ \ ll \-2 . 2 \ -3/4+1/100 

(4.6) 

||e^V(t2)|L2 <C{iEn)-' + ^on's)-\ 

loc 

These estimates are used to estimate the local contributions of e^'-'riE{t2) in (3.19) 
(see the following remark) and in the proofs for the and L^^^ estimates for r]E{t) 
in Lemmas 5.2 and 5.3. of [P!^ . 
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Remark: In [|r^, is used to estimate fi^\t) = e~*'^^*^*^^r/(t2), with ^7(^2) = 

e*©(t2)f/^^(^^2), and A being a self-adjoint perturbation of —A + V, C = U^^{—iA)U, 
U a bounded operator. Since e'^^ = U~^e~^'^^U, 

Thus we can choose freely to estimate either A or C Since e^^~^'-^^^r]E{t2) is easier to 
estimate than e"*'^^*"*^'' 7/(^2) by using the Duhamel's expansion, we state all conditions 
in terms of C. 

Because we only need ([4.6|), the same proof in actually gives the following 
stronger result. 

Theorem 4.1 Suppose that tp{t2) is close to a nonlinear ground state Qso^*®" in -^foc" 
norm, and suppose that in the decomposition ( [4.3|) with E = Eq one has WQeoWl'^ = 
nt.,r^n< no, \\Ceo\\ < ^o^, \a\ + Hr/eJI^a < eln^ . 

loc 

If for all E close to Eq with \E — Eq\ < ej^n^ the dispersive part riE{t2) in the 
decomposition (|4.3|) satisfies the estimates (|4.6|) , then the conclusion and the proof of 
Theorem 1 in / [7^/ remain valid. In particular, there is a frequency E^o with \Eoo — 
Eo\ < ^0^^ '^'^^ ^ function Qit) = —Eoot + 0(log(t)) for t E [t2, 00) such that 

\\m - Qi^ooe^^'^^^L. < C2 {{en)-' + jon\t - t2))"'^' 

loc 



for some constant €2- 



Suppose that the initial data ip{t2) is decomposed as in (|4.3| ) with the frequency E 
chosen so that the coefficient a vanishes. Suppose that, in addition to the previous 
assumption that ( [4.6|) holds for all frequency E with \E — Eo\ < s^n? , the excited state 
component ( satisfies ( [4.4|) - ([4.5|) . Then the lower bound 

Ci {{en)-' + ^on\t-t2)y'^' < - Qi^^e^^^ | 

holds as well. 



loc 



The merit of this modification is that we do not need the initial data to be local- 
ized. We only need to know that its dispersive part is "outgoing" in certain sense. 
Notice that the condition on the size of the excited component C is a simple conse- 
quence of the estimates (p.6|), (|3l8| ) and ( p.4| ). To see this, we first pretend that the 
size of ( is given by 7/(^2) and the size of the ground state component is given by x(t2)- 
Then the condition ( |4.4| ) is just a simple consequence of ( p.8| ). Since the difference 
between the decompositions ( |1.12| ) and ( ^4.3| ) are higher order terms, the condition 



(O) is easy to check. Therefore, Theorem p..l| follows from the following Lemma: 

34 



Lemma 4.2 Let %l){t) he the solution of ( |1 . 1| ) in Theorem |i.i| and t2 he the time in 
Proposition \3. 1\ . Let Eq = E{t2) he the unique energy such that in the decomposition 
( [4. 3D the coefficient a vanishes. Then for all E close to Eq, i.e., \E — Eq\ < Ce^n"^ , 
we have for all t > t2 



e 



loc 



Notice that, since t2 ~ by Proposition we have (£:n)^^ + 7on^(t — ~ 

n?t for all t >t2, no matter t > 2t2 or t < 2t2- Hence (|4.71 ) implies ( [4 .61) with a big 
margin 

Proof of Lemma |4.2| 

We have the two decompositions at t = ^2 

ij{t) = x{t)^o+QMt))+m 

(4.8) 

= [Qe + aE{t)RE + CE{t) + riEit)] e*®^W 
Since E will be fixed for the rest of this proof, we shall drop all subscripts E. Hence 

at) + r/(t) = [x{t)(t)o e-*®(*) - Qt] + [Qi{y{t)) + m] e"*®^*^ - aE{t)RE . (4.9) 
Thus we have 

r^{t2) = P/ {[x(t2)0oe-^®(*^) - Qt] + [^1(^/(^2)) + ^(^2)] e'^^^*^)} = 7^0,1 + ^/0,2 , 
where 

r/0,1 = P/ {e(t2)e-'®(*^)} , 

r/0,2 = P/ {(x(t2)0oe-*®(*^) -Qi,) +Qi(y(t2))e-^®(*^)} . 

Note that 770,2 is a local if^ function and is bounded by 0(n'^), i.e., ||77o,2||y < Cn^. 
Therefore we have 

||e(*-*^)S,2|L. < Cn' {t - t2)-"' , ||e(*-*^)S,2|L. < C7n^ (t - t2)-"' . 

loc 

(4.10) 

Hence e^^T]Q 2 satisfies the desired estimates with a big margin. 

We now focus on the non-local term 770,1 = {^(^2)6"*®^*^-'}. Recall ^(^2) is 
bounded by 2n in (|3.2|) and by 4(72^^"^/^ in L^^^^ from Proposition |3.1| . In particular, 
we have ||77o,i|Il2 < C 11^(^2)11^2 < Cn. Thus, by Lemmas 2.6 and 2.9 of [0, we have 
||e(*-*2)^77o,i||^2 < Cn. 
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For convenience of notation, we write 

L = -iHo + iVi + iV2C , 

where Vi = 2AQ|;, V2 = and C denotes the conjugation operator. By Duhamel's 
principle, 



We now substitute (|2.14| ) with t = t2, i.e., 

^(^2) = e-**^^oeo + /*' e-*(*^-")^« Pc^« G^(r) rfr, 
Jo 



(4.11) 



into the above equation. We have 



e^'-''^%,i = E, + E2 + E^ + E^ 
^2= / e(*-^)^P/z(ri + \/2C)e-'^^°e-'®(*2)eoc^s 

Jt2 

Jt2 Jo 



where 



The only estimates we need here from §3 are: (cf.( p.l4| )) 



\My<An, \\Gi:{r)\\^,<Cn', (0 < r < ^2), ^2 ~ V 



(4.12) 



From the linear estimate Lemma 2.2 we can estimate Qi by 



For $72; since ^2 > 1? we have 



11^21^4 <C f\t- st""' lleolLs/r ds < Ct-'/' Uohs/r 

Jt2 
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If t > ^2 + 1; we bound its L^^^-norm by for s E [t2,t-l] and by L'^ for s E [t-l,t]. 
Thus we have 

<cr^/«lieolL8/r . 

If if: < )f:2 + 1; we use only L'^ norm for the whole interval [^2, t] and the same estimate 
holds. 

For il^, notice that ^l^ = Pj^ e~^^^~^'^^^'^J where J is the integral in ( [4.11|) , 

J= e-^(*2-")^° Pc^° G^(r) dr = ^(^2) - e-^^^-^^^o • 
Jo 

Since it is the difference of two functions of order n, \\J\\i^2 < Cn and ||r23(t)||j^2 < 
Cn for all t. 

Suppose that t > 2t2. Since HG^H^i < Cn^ by ( p. 141 ), we have 



Jo io 

Interpolating, we have 

Since L\^^-mixva. is bounded by L°°-norm, we conclude that 

11^^311^2 < CnH^t-^l'' < Cn't'Jh-^l'' . 



Suppose now that t < 2^2- From similar arguments we have 

11^311^4 <C r \t - T\~^/^n^ dr < C r \t2-T\-^'^n^dT = Ct 
Jo Jo 



1/4 3 

2 n 



-12-1 r^2 

It2-1 

"12-1 rt2 



< 



C [' 1^2 - T\-^^^nUT + C \t2- r|-3/Vdr < Cn^ + Cn^ 

Jo Jt2-l 



We now estimate ^4. Since HG^H^i < Cn^ by ( p.l4| ), we have 

"t rt2 



IWIl^ <C [ \t- s\-^/^ \s- t\-''/^ drds 

Jt2 Jo 

<C f \t-s\-^'%s-^/^n^ds 

Jt2 

<Ct2n^ f \t-s\-^/h-^/Hf"'ds 

Jt2 

<CnHT f\t-s\~^"{s)-'" ds 
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We can bound the Li^^-norm by for s G [t2, ^ — 1], and by for s E [t — l,t] 
(if t < ^2 + Ij we use only norm for the whole interval [t2)^])- Thus we have the 
bound 

\W\lI^< {c j'~\t-s\-'^/'' + C j' i^f^' \s-r\-^'^\\G^{T)\\^, dr^ ds 

<C f \t- s\-'"\t2S-''^^rv')ds + C f \t- s\-^l\t2S-^'^rv')ds (4.13) 



Jt2 Jt-l 

The second integral in ( [4.13|) is bounded by Cn^t2t^^^'^. The first integral, when 



t >2t2, is bounded by 

■>t/2 nt-1 



< Cn% [ \t- s\-'/'s-'/' ds + Cn% [ \t - s\-'/'s~^/^ ds (4.14) 

Jt2 Jt/2 

<CnH2r'/' f'\-^'^ds + CnH2t-^'^ f \t-s\-'^/^ds 

Jto Jt/2 



't2 Jt/2 

< CnH2t-^/%'^' + CnH2t-^/^ 



On the other hand, ii t2 + 1 < t < 2^2, then the first integral in ([4.13 ) is bounded by 



t-i 



CnH2 / \t-s\-^lh-^'^ds 

Jt/2 



which is the second integral in ( |4.14| ) and is bounded by Cn^t2t Thus, for all 

te[t2,T], 

11^^411,2 <CnHl^h-^/^ + CnH2t-^/' 

loc 

Summarizing, we have 



loc ' ' loc I. ) 

Since ^2 ^ e^^n""^, equation ( |4.7| ) holds if n is sufficiently small. We have proved 
Lemma 
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